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ABSTRACT 

This  report  presents  certain  relations  among  the  completely 
integrable  Hamiltonian  systems  introduced  by  Gel' fand  and  Dikii. 

These  relations  generalize  a formula  of  A.  Lenard  linking  the 
higher-order  Korteweg-deVries  equations,  of  which  the  Gel ' fahd-Dikii 
Systems  are  a generalization.  The  general  form  of  the  relations, 
which  connect  the  various  isospectral  deformations  of  linear  differ- 
erential  operators,  is  described,  and  two  examples  are  given 
explicitly. 
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SIGNIFICANCE  AND  EXPLANATION 

In  recent  years  several  nonlinear  partial  differential  equations  of 
applied  mathematics  have  been  discovered  to  have  the  extraordinary  property 
known  as  complete  integrability : that  is,  roughly  speaking,  they  possess 
the  maximum  number  of  constants  of  motion  possible  for  the  type  of  system 
considered.  These  equations  arise  in  the  study  of  shallow  water  waves, 
accoustic  waves  in  plasmas  (Korteweg-deVries  equation)  nonlinear  optics, 
Josephson  junction  theory  (sine-Gordon  equation) , other  plasma  phenomena 
(nonlinear  Schrttdinger  equation) , and  other  areas.  The  complete  integra- 
bility property  - again,  roughly  speaking  - allows  unusually  explicit 
solution  of  these  equations. 

These  equations  are  undoubtedly  very  special.  It  seems  important  to 
understand  the  place  occupied  by  these  special  systems  in  the  general  ag- 
gregate of  partial  differential  equations.  Recently,  I.  M.  Gel' f and  and 
L.  A.  Dikii  have  succeeded  in  isolating  several  key  facets  of  the  structure 
of  these  systems,  simultaneously  constructing  families  of  hitherto  unknown 
systems  with  the  complete  integrability  property.  This  construction  reveals 
slightly  more  clearly  the  nature  of  these  systems. 

Each  of  the  aforementioned  partial  differential  equations  - Korteweg- 
deVries,  sine-Gordon,  nonlinear  Schrodinger  - is  embedded  in  a heirachy  of 
equations,  each  having  the  complete  integrability  property,  and  each  related 
to  the  others  by  a common  rule,  discovered  in  the  case  of  the  Korteweg-deVries 
heirarchy  by  A.  Lenard.  This  relation  allows  a transparent  derivation  of  some 
of  the  extraordinary  structure  of  the  Korteweg-deVries  equation.  This  report 
presents  a similar  rule  relating  the  various  hierarchies  of  systems  intro- 
duced by  Gel 'f and  and  Dikii. 
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W . Symes 

§0.  Introduction 

In  this  short  note  we  prove  the  existence  of  certain  relationships  amongst  the 
Hamiltonian  systems  constructed  by  Gel'fand  and  Dikii  in  (1).  These  relations  generalize  a 
formula  of  A.  Lenard  linking  the  higher-order  Korteweg-deVries  equations,  of  which  the  sys- 
tems in  [1]  are  a generalization.  We  refer  the  reader  to  [2] , 13]  for  a description  of 
Lenard' s result,  which  is  also  reproduced  as  an  example  at  the  end  of  this  paper,  and  for 
an  explanation  of  its  importance  in  the  theory  of  the  Korteweg-deVries  equation  and  its 
higher-order  relatives.  We  anticipate  similar  applications  to  the  Gel ' fand-Dikii  systems. 

Mark  Adler  has  derived  the  generalized  Lenard  relations  of  this  paper  in  somewhat 
different  ways,  first  for  hierarchies  of  systems  related  to  the  Boussinesq  equation  and  to  a 
certain  fourth-order  eigenvalue  problem,  [3],  then  in  general  for  the  Gel'fand  Dikii  systems 
[4].  His  work  will  appear  elsewhere. 

I would  like  to  thank  Mark  Adler,  Charles  Conley,  and  Neil  Fenichel  for  a number  of 
illuminating  conversations  on  this  subject. 

§1 . Formal  Isospectral  Deformations. 

Let  be  a linear  ordinary  differential  expression  of  order  n : 

L=  Dn  + l q Dk  , D = - i ± . (1) 

k=0  K dx 

The  domain  of  definition  of  the  coefficients  q^,  which  are  supposed  infinitely 

differentiable,  is  some  open  interval,  u , finite  or  not,  in  the  real  line. 

A formal  isospectral  deformation  of  L is  a specification 
n-2 

l = l q.  D = IPrL]  (2) 

k=0 

where  P is  a differential  expression  whose  coefficients  depend  polynomially  on  q^,...,^  ^ 
and  their  derivatives,  having  the  property  that  the  commutator  appearing  on  the  right-hand- 
side  of  (2)  is  of  order  n-2  or  less.  Thus  (2)  may  be  regarded  as  a collection  of  ( n— 1 ) 
partial  differential  equations  for  the  coefficients  q^,  where  the  dot  is  interpreted  to  mean 
differentiation  with  respect  to  a ("time")  parameter.  Peter  Lax  has  shown  that,  if  L is 

provided  with  a suitable  function  space  domain,  becoming  an  operator  and  attaining  a spectrum, 
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then  a solution  of  such  a system  of  partial  differential  equations  represents  a deforma- 
tion of  L preserving  its  spectrur.  (See  [5]). 

For  exairple,  let  n = 2.  Then 

L =D2  + qQ  , P = D3  + ^ qD  + | (Dq) 

and  the  equation  (2)  is  equivalent  to  the  Korteweg-deVries  equation  for  the  coefficient  qQ: 

_i  ^0  = 4 °3q0  + I q Dq0- 

Of  course,  the  existence  of  a suitable  expression  P is  clearly  a local  matter,  having 
nothing  whatever  to  do  with  selection  of  a domain  for  L.  Gel’fand  and  Dikii  present  a con- 
struction in  tl]  for  all  such  P , based  on  the  local  algebra  of  symbols,  which  is  the  formal 
side  of  the  calculus  of  pseudodifferential  operators.  We  give  a very  brief  review  of  this 
construction  in  §12,3,  and  4.  All  of  the  statements  in  §§2,3  are  formal  counterparts  of 
results  of  Seeley,  [6],  and  proofs  are  also  sketched  in  [1],  We  therefore  omit  the  proofs. 

§2 . Symbol  Algebra 

By  intentional  abuse  of  notation,  denote  also  by  L the  symbol  of  L , 
n-2  . 

L[q,U  = C"  + l qv  C (3) 

k=0 

where  we  have  emphasized  the  dependence  of  the  symbol  L on  the  coefficient  vector 
q = (q0 qn-2>  * 

More  generally,  a symbol  is  a formal  sum 

00 

A(x,£)  = l A£(x,C)  (4) 

1=0 

where  A^  is  smooth  and  complex-valued  in  U * ®\{0},  and  homogeneous  of  degree  d^  in  6 

for  large  | £ | , with 

d„  > d,  > ...  > d„  > 

0 1 l 

a sequence  of  real  numbers  tending  to  -■»  . 

Symbols  are  considered  the  same  if  they  agree,  term-by-term,  with  each  other  for  5 out- 
side some  sompact  neighborhood  of  0 e K . In  that  case  we  will  write  the  sign  of  equality, 
remembering  that  it  means  "for  large  | C I " - 

The  symbols  form  a module  over  C°°(U)  , with  the  module  operations  defined  pointwise. 

They  form  an  algebra  over  I,  with  the  multiplication  operation  o defined  by 


(5) 


r 1 ^ v v 

(A,B)  H>AoB=  l — , 3 A D B 

v>0 

3 3 

where  3 = — and  D « -i  - — are  applied  term-by-term.  Note  that  the  sum  on  the  right  of 
(5)  is  not  in  the  canonical  form  (4) ; however,  only  finitely  many  products  of  the  same  homo- 
geneous degree  appear  in  (5) , and  rearrangement  into  the  form  (4)  is  easy. 

Note  that  the  class  of  symbols  all  of  whose  homogeneous  pieces  have  integral  degree, 
and  the  class  of  symbols  which  for  large  | C I can  be  written 

A = I \ tq.si 

i>0 

de 

with  A,  [q,51  = a^  [q]  (,  and  a^  [q]  a polynomial  in 
a polynomial  in  • • • ,<Jn_2  and  x-derivatives)  , 

subalgebra  will  be  called  the  class  of  q-symbols. 

The  order  ord  A of  a symbol  A is  the  homogeneous 
geneous  degree  appearing  in  an  expansion  (4)  - that  is. 

Note  that 

ord  (A  o B)  <_  ord  A + ord  B 

and 

ord  [A,B1  _<  ord  A + ord  B - 1 

where 

[A,B]  = A«B-BoA. 

Note  that  differential  expressions  such  as  (1),  and  more  generally  differential  expressions 
whose  coefficients  depend  polynomially  on  q and  its  x-derivatives,  correspond  under  the 
replacement  D 5 to  polynomial  (in  O q-symbols.  Moreover,  this  correspondence  commutes 
with  the  formation  of  products  and  sums.  Therefore,  Lax  equations  (2)  are  in  1-1  correspon- 
dence with  polynomial  q-symbols  P for  which  ord  [L,P1  < n-2. 


q and  its  x-derivatives  (that  is, 
both  form  subalgebras.  The  latter 

degree  of  the  term  of  highest  homo- 
ord  A = dQ  in  the  notation  of  (4) 
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S3.  Complex  Powers 


Define  the  resolvent  symbol  R(A)  for  L by  the  equation 
R(A)  ° (L-A)  - 1 . 

According  to  the  definition  of  the  product  (5),  we  can  re-write  this  equation  in  the  form 

n-2 

<Sn-A)  R ( A ) = 1 - R(A)  • ( l q.  5 ) • 

k=0  K 

From  this  formula  it  is  easy  to  see  that  R(A)  admits  an  expansion  in  the  form 

R(A)  ■ [ R (A.O  (6) 

no- 
where R^  is  homogeneous  of  degree  -£  in  £,  /A  , and  is  defined  recursively  by  the 

formulae 


R (A)  = (5n-A)-1 
n 


Rn+1(A)  , 0 . 


i-n-2  n-2 

*,<*>“-  l l 


j=0  k=max(0,2n+j-l) 


ek(C"  -A  )-1  -2n-j+£+k  -2n- j+£+k 

(-2n- j+£+k) 1 n+j  qk 


In  particular, 

£-n 

R,(X)  = l B (c"  + A)-"  , £>  n+2  (7) 

* m=2 

where  B.  is  a homogeneous  polynomial  in  5 of  degree  nm-£  whose  coefficients  are  poly- 
£m 

nomials  in  q emd  its  x-derivativesr. 

Observe  that,  by  virtue  of  (7),  the  expansion  (6)  can  be  rearranged  for  large  |?| 
into  the  form  (4),  showing  that  R(A)  is  a q-symbol.  R(A)  should  be  regarded  as  a local 
version  of  the  resolvent  of  L . Indeed,  R obeys  the  resolvent  equation 
R(X)  - R(w)  = (A  - v)  R(A)  o R(U). 

In  further  analogy  to  the  usual  theory  of  the  resolvent,  we  use  R ( A)  to  define  complex 
powers  of  the  symbol  L . 

For  the  remainder  of  this  paper,  define 
A 3 = exp(s  log  A),  A,s  e <1 

where  the  principal  branch  of  log  is  selected,  with  the  branch  cut  down  the  negative  imaginary 
axis. 

Define  the  symbol  L , for  complex  s , by 

LS  = ~r  f d>  AS  R(X)  *(t)  (8) 

The  contour  goes  up  the  ray  Re  A=  1/2,  around  the  semicircle  | A | = 1/2  counter-clockwise, 
and  back  down  the  ray  ReA  = -1/2.  The  function  x<C)  is  smooth  and  satisfies 


for  R < A ) in  (8) 


x (O 


r i . id  1 1 

V 0 . k|  < 1/ 


;|  < 1/2 

The  integral  (8)  is  evaluated  by  inserting  the  expansion  (6) , (7) 
and  integrating  term-by-term,  using  the  Residue  Theorem.  Each  integral  converges  and  admits 
evaluation  by  residues  for  Re  s < -1:  We  obtain 

(9) 


where  ord  A 


LS=  I A (s) 

p=0  y 

n Re  s - p , 


A (s)  = (C  >' 

P 


<C‘V"“  f B (-1) 

m=2  P+n'm 


m-1 


U 


(10) 


with 


m-2 

l (s-j) 
j=0 


f s ) = ± 

Tit-l'  (m-1)  ! 

According  to  the  convention,  that  symbols  are  identified  when  they  agree  for  large  | £ | 


the  above  results  (9) , (10)  are  independent  of  the  choice  of  x » which  appears  in  (8)  merely 
to  rule  out  poles  on  the  integration  contour.  All  mention  of  x in  the  following  will 
therefore  be  suppressed. 

g 

Thus  L is  a q-symbol,  of  order  n Re  s,  for  Re  s < -1.  Formulae  (9) , (10)  allow 
• • s 

continuation  of  L as  an  entire  function  of  s . Himicing  arguments  in  [6) , one  easily 
shows  that  the  Resolvent  Formula  implies  the  semigroup  properties 


L = 1 , L 
for  any  complex  s,t, 


L , LS  Lfc 


P s+t 


[LS  , LfcJ  = 0 . 


In  particular 


§4 . Lax  Pairs 

Let  m be  any  positive  integer.  Then  Lm^n  is  a q-symbol  of  order  m , whose  homo- 
geneous pieces  have  integral  degree.  Set 
Lm/n  = Pm  , N1" 


where 


= [ A ( — ) 

c p n 


r . in  . 

y a < - > . 

n n 


p=m+l 


p n 


Thus  Pm  is  a polynomial  q-symbol,  and  ord  Nm  -1  . If  m is  not  divisible  by  n , it 

is  easy  to  check  that  A ( — ) ^ 0 for  p > irH-1  . 

p n — 

, m/n , 

Because  [L,  L ] = 0 , 

(Pm,  L]  = [L,  N1"]  . 

Since  the  left-hand  side  of  this  equation  is  a polynomial  q-symbol,  so  is  the  right.  On  the 
other  hand;  since  the  right-hand  side  has  order  <_  n-2,  so  does  the  left-hand  side.  Thus 
£ = [pm,  LJ 

is  a Lax  equation  for  each  positive  integer  m . This  is  Gel'fand  and  Dikii's  construction 
of  Lax  pairs. 

A simple  degree-counting  argument  shows  that  P™  is  essentially  unique  polynomial 
q-symbol  of  order  m , whose  commutator  with  L is  of  order  n-2. 


§5.  Generalized  Lenard  Relations 

Denote  by  Xm  the  polynomial  q-symbol 
Xm  s [Pm,  L]  = [L,  Nm]  . 


objects  of  this  section  is  to  relate  X to  X 
p"*"  = (L  . Lm/n)  . 


xm 

, th 

the  m 

Lax  vector 

field  for 

L . The 

m+n 

. Since 

m/n  + 1 

L — 

L o Lm/n  , 

we  have 

= L o P1"  + (L  ° n"1) 


where  the  subscript  "+"  signifies  the  sum  of  parts  of  positive  homogeneous  degree,  i.e. 


the  polynomial  part.  Thus 
m+n 


L o X1"  - [L,  (L  o Nm)+]  . 


(11) 
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The  equation 


n m+k-i-1 

XT  = I l [A ( qt  tV 


k=I+2  p=m+l 


p n ^k  ” ‘l 


(where  the  subscript  i denotes  the  part  of  homogeneous  degree  t)  shows  that  xm  depends 

only  on  A^(  — ),  m+1  <_  p <_  m+n-1.  We  shall  show  that  Equation  (11)  can  also  be  rewritten 
nw-n  m 

to  express  X in  terms  of  Ap(  - ) , m+l^p<  m+n-1.  This  will  be  our  generalized 


Lenard  relation. 


In  fact. 


(L  o Nm)  '=  A ( S.  ) + ... 
+ m+n  n 


where  the  terms  denoted  by  dots  involve  only  Ap(2)ftn+1^p^  rn+n_i.  0n  the  other  hand 
according  to  the  result  of  §4, 

0 = [L,  Nm]_1 

= I ( " ) Dn+1‘P  a < 5 , 

n+l-p  m+p  n 


l l < v+L>  fa. 


k=0  p=l 


k+l-p'  ^k 


, m , m _ k+l-p  , 

am.  ( ~ ) - z ( ~ ) (-D)  r q,  } 
m+p  n m+p  n k 


where  we  have  written  A ( ^ ) = a (51)  This  can  **  re-written 

p n p n 


-n  D a ( — , = £ ( " ) Dn+1_P  a ( * , 

m+n  n L.  n+l-p  m+p  n 


n-2  k 


+ I I < * > (qk  Dk+1-P  a ( ^ ) - a ( ” ) (-D)k+1-p  ql  (13) 

, r,  „=n  k+i-p  k m+p  n m+p  n k 


k=0  p=l 


The  r.h.s  of  (13)  is  an  exact  derivative,  so  we  can  express  a ( — ) in  terms  of 

m+n  n 

am+p  1 n ^ ^ — p — n-1-  0,1  the  other  hand,  we  don’t  even  need  to  do  that,  since  only 
DAm+n  ^ n * appears  in  [L , (L  o Nn)  + J.  In  any  case,  we  have  expressed  xm+n  in  terms 
of  A_  ( — ) > 1 < p < n-1,  as  was  desired. 


§5.  The  Hamiltonian  Formalism 

In  order  to  write  the  relations  derived  in  the  last  section  in  the  form  in  which  the 
Lenard  relations  are  usually  presented,  we  introduce  the  Hamiltonian  formalism  of  11]- 

Suppose  that  F is  a polynomial  in  q and  its  derivatives.  Define  the  formal  varia- 
tional derivative  of  F by 


£-=  I <-D>; 


3F 


j=0 


3 (D^) 


3F  . 3_F  3 F 

3<30 3qn-2 


) 


. m . n 

a i — ) = — — 

p n m+n 


S 


In  [1]  it  is  proven  (Equation  (21))  that 
p-m-1 

i ‘pti>‘-d>v  , 

v=0  q?-m-v-l 

for  p = m+1 , . . . , nH-n-1 . 


m+n+1 


( ^ ) 
n 


(14) 


Using  (12)  and  (14)  , it  is  easy  to  express  the  Lax  vector  field  X in  the  form 


(X™ Xm  ,)t  = -7-  J H 

0 n-2  m+n  6q  m 


(15) 


where  = am+n+^  ( ~~  ) » and  J is  a certain  matrix  of  differential  expressions,  whose 
coefficients  depend  polynomially  on  q and  its  derivatives. 

In  [1]  it  is  shown  that  J defines  a synplectic  structure  on  the  space  of  coefficient 
vectors  q , in  a certain  local  sense,  and  that  the  Lax  vector  field  is  a Hamiltonian  vector 
field,  by  virtue  of  (15) . 

Now  Equations  (11),  (12),  (13),  (14)  and  (15),  allow  us  to  write  the  relations  derived 
in  the  last  section  in  the  succinct  form 


J f H_  = -£-K  f H 

m+2n  oq  m+n  n+m  oq  m 


(16) 


where  K is  another  (n-l)x(n-l)  matrix  of  differential  expressions  whose  coefficients 
depend  polynomially  on  q and  its  derivatives.  Equation  (16)  is,  for  n=2,  the  form  in 
which  the  Lenard  relations  are  usually  stated.  The  general  form  of  J is  very  comp-icated, 
and  that  of  K even  more  so.  Rather  than  exhibit  these  general  forms,  we  compute  J and 
K explicitly  for  n=2,3  in  the  next  section. 
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§ 7 . Examples 

2 

I.  n = 2 : L = F + q^  , and 
0 


m , m . 4 6 . m , . 

X = 2 D a ( - ) = —z  D t — a _ ( - + 1) 
0 m+1  2 m+2  m+3  2 


so  in  this  case  J is  the  1X1  matrix  2D. 

Equation  (11)  yields 

m+2  2m  m ^ 3 . m ( 2 , m v . m . 

x0  = D X0  + q x0  - 20  am+l  < 2 > - D am+2  ‘ 2 > + Vl  ( 2 > °q 


Equation  (13)  yields 


Da  „ = -1/2  Da,. 
m+2  m+1 


X^2  = ( ^D3  + 2 qD  + (Dq)>  am+1  ( f ) 


So  K is  the  lxl  matrix 

K = ;D3  t qD  t D 

2 ^ Oq 

This  is  the  expression  discovered  by  Lenard. 

II.  n = 3.  Set  H = a ( ^ + 1)  . 

m m+3  m+4  3 


m+1  6q 


am+2  6q  ^ 6q 


according  to  formula  (14)  . 


The  Lax  vector  field  is 


Xl"  = 3«Dam+l  < I > + 3(d2  am+l  < ! > + D a2  < f » 


From  (13)  we  obtain 

D V3  < ! > = - I (3C2  am+2  < I > + D'  am+l  < ! > + (D  ‘ V a^  ( 


and,  after  tedious  computation 


w|3 


where 


KQ0  = - | 1)5  * 3 1)3  " ‘ j'^l  1)3  ' 3 qilD  ° + D ° + <Dqo,D 

Koi  = °4  - 2qi  °2  + 3qoD  ' Dqo 

K1q  = -D4  + 2 D2  o qx  - 2 ( DqQ ) + 3qQD 
KU  = 21)3  + qxD  + D o qi  . 
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